Nernst effect anisotropy as a sensitive probe of Fermi surface distortions 

from electron-nematic order 
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We analyze the thermoelectric response in layered metals with spontaneously broken rotation 
symmetry. We identify the anisotropy of the quasiparticle Nernst signal as an extremely sensitive 
probe of Fermi surface distortions characteristic of the ordered state. This is due to a subtle interplay 
of different transport anisotropics which become additionally enhanced near van-Hove singularities. 
Applied to recent experiments, our results reinforce the proposal that the underdoped cuprate 
superconductor YBa2Cu306+i displays such "electron-nematic" order in the pseudogap regime. 



Spontaneous breaking of lattice rotation symmetry 
due to electronic correlations is currently in the focus 
of intense interest, most prominently in cuprate high- 
temperature superconductors such as YBa2Cu306+(5 
(Refs. 1-3) and in the metamagnetic metal Sr3Ru2 07 
(Ref. 4). In analogy to liquid crystals, a phase with bro- 
ken rotation (but preserved translation) symmetry has 
been dubbed electron nematic.^ 

In cuprates, electron-nematic order has been discussed 
early on as intermediate phase which occurs upon melt- 
ing of a uni-directional charge-density- wave ("stripe") 
phase. Microscopically, it is one of the known insta- 
bilities of the two-dimensional (2d) Hubbard model. ^"^"^ 
The first clear-cut signature of electron-nematic order 
in cuprates was found in neutron scattering experiments 
on YBa2Cu3 06.45,^ where the spin-fluctuation spectrum 
was found to develop a distinct anisotropic incommen- 
surability below a temperature T of about 150 K. Ear- 
lier transport measurements on YBa2Cu306+5 (Ref. 11) 
also detected resistivity anisotropics Pa/ Pb of order 2 in 
the underdoped regime, but remained less conclusive. 
The reason is that the crystal structure of YBa2Cu306+5 
contains CuO chains which break the otherwise tetrag- 
onal symmetry of the Cu02 planes. In an order- 
parameter language, this implies a small field which cou- 
ples linearly to the nematic order parameter. This has 
two main effects: (i) potentially existing nematic or- 
der will be aligned and (ii) a nematic ordering tran- 
sition will be smeared out. While (i) enables observ- 
ables to show a macroscopic anisotropy, which might 
otherwise be masked by domain formation, (ii) implies 
that electronic nematic order and purely structural ef- 
fects cannot be sharply distinguished. Remarkably, lo- 
cally broken rotation symmetry has been found on the 
surface of Bi2Sr2CaCu208-i-5 and Ca2-a;Na2,Cu02Cl2 us- 
ing scanning tunneling microscopy.^^ In Sr3Ru207, ne- 
matic order is a candidate explanation for the low-T 
phase which masks the metamagnetic critical endpoint 
at around 8T.^'^'^ Resistivity anisotropics have been de- 
tected here, but a full picture has not yet emerged, be- 
cause the rather unusual thermodynamics near the low-T 
phase boundaries is not understood. 

Very recent measurements^'^ of the Nernst effect, 
that is the transverse voltage induced by a thermal gradi- 



ent in the presence of a magnetic field, have uncovered a 
surprisingly large anisotropy in YBa2Cu306+5 below the 
pseudogap temperature. Consequently, this has been in- 
terpreted as evidence for broken rotation symmetry, but 
a theory which links this anisotropy to a particular or- 
dered state was not available. Previously, the Nernst 
effect had been extensively studied as a probe of cuprate 
pseudogap physics, but mainly interpreted as sign of 
precursor superconducting fluctuations. Only recently, a 
second contribution to the Nernst signal has been iden- 
tified and associated to quasiparticle physics. Theoret- 
ically, the quasiparticle Nernst effect in the presence of 
antiferromagnetic^^ and stripe^^ order has been investi- 
gated. 

In this paper, we analyze, for the first time, the spa- 
tial anisotropy of the Nernst effect. We show that an 
enhanced and strongly anisotropic quasiparticle Nernst 
signal arises from the Fermi surface distortions which 
accompany broken rotation symmetry in metals. Sur- 
prisingly, small anisotropics in the kinetic energy can be 
easily enhanced by an order of magnitude in the Nernst 
signal, due the interplay of the anisotropics of conduc- 
tivity and thermopower. A similarly large anisotropy 
does not arise from stripe phases. Hence, we establish 
the Nernst effect as a very sensitive and unique probe of 
electron-nematic order. Applied to the measurements 
on underdoped YBa2Cu306+5 (Refs. 2,3), we confirm the 
presence of the long-sought electron-nematic phase in the 
pseudogap regime. 

Quasiparticle model. We shall calculate the thermo- 
electric response in a one-band model of quasiparticles. 
For concreteness, we work on a two-dimensional (2d) 
square lattice with a tight-binding dispersion 
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Such a 2d model is appropriate not only for cuprate su- 
perconductors, but (with some modifications) for a va- 
riety of other layered materials with tetragonal lattice 
structure, including Sr3Ru207. To simplify calculations, 
we shall neglect inter-layer coupling, but will comment 
on its effect towards the end. 

Nematic order distorts the band structure, which, on 
the mean-field level, can be captured by anisotropic hop- 
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ping parameters.*'^" Here we shall focus on the case of 
a d^2_y2 electron nematic, as may arise from a Pomer- 
anchuk instability in the I — 2 channel or as precur- 
sor of a stripe phase. We introduce an anisotropy pa- 
rameter e, such that the hopping matrix elements obey 
tia;,y = (1 ± e/2)ti and ta^^y = (1 ±e/2)t3. For e ^ 0, the 
lattice symmetry is thus broken from C4 down to C2. 

Transport theory. The linear thermoelectric response is 
captured by three conductivity tensors a, a, and k, which 
relate charge current J and heat current Q to electric 
field, E and thermal gradient, VT vectors: 

(;J>(/.t)(-i)^ (^) 

The electrical field induced by a thermal gradient in the 

absence of an electrical current can be expressed by the 
linear response relation E = — i?VT, and Eq. (2) together 
with J = yields E = a-'^aVT. Therefore, the Nernst 
signal {)yx, defined as the transverse voltage Ey generated 
by a thermal gradient V^T, reads 

q ^XX^yX ^yX^XX /ry, 

^xx^yy O^xy^yx 

and "dxy is obtained from x ^ y. For a magnetic field 
B = Bz in z direction, the Nernst coefficient is usually 
defined as Vyx = 'dyx/B, which tends to become field- 
independent at small B. We employ a sign convention 
such that the vortex Nernst coefficient is always positive 
(formally Vxy = —"^xy/B, Ref. 20). In general, the Nernst 
coefficient can be negative or positive, for example if it 
is caused by the flow of charged quasiparticles. 

The simplest description of quasiparticle transport is 
via the Boltzmann equation in the relaxation-time ap- 
proximation. We assume a momentum and energy- 
independent relaxation time tq, as is appropriate for elas- 
tic impurity scattering at low temperatures. Then, the 
electrical and thermoelectrical conductivities read: 

- The ■^de{\.f^^>'''^V'^dkx ''SfcJ' 

<rxx - -2ero;^^K) , 

a - -2^V^.4.^^-.^^1 (4) 
''^^ - ^ he ^de{\,f''\^'^dkx "^dkyY 

where e = |e| is the electron charge, Vk the velocity, 
and the Fermi function. In this approximation, the 
Nernst signal is linear in both T (at low T) and tq. In 
the following, we shall be interested in the magnitude 
and anisotropy of the low-T Nernst coefficient, i.e., we 
shall use Eqs. (4) to calculate Vxy/T and Uyx/T. 

Analytical considerations. On general grounds, the 
Hall conductivities obey Gxy = —Oyx independent of 



the crystal symmetry. Such a relation does not hold 
for ctxy^yx in general, however, in the low-T limit the 
Mott relation can be derived from Boltzmann theory,^^ 
d oc da/dn, implying that axy = —(Xyx- Any anisotropy 
in the Nernst signal is therefore arising from axx—o'yy ^ 
and OLxx ~ Oiyy 7^ - both will have a piece which is lin- 
ear in e. In general, the Nernst anisotropy at low T will 
depend on the anisotropies of the Fermi surface, of the 
Fermi velocities, and of the scattering rate (the latter 
anisotropy is not taken into account in our calculation). 

It is worth noting one limiting result: For the sim- 
plest anisotropic dispersion, ek = {k'^/iT^x + ky/my)/2 
with niy^x = (1 i e/2)m, the anisotropies of a and a 
obtained in the Boltzmann framework (4) cancel in the 
Nernst signal, i.e., the Nernst signal remains zero due 
to Sondheimcr cancellation.^'' (t^ = is an artifact of an 
energy-independent tq.) However, we shall obtain a large 
anisotropic Nernst response for non-parabolic dispersion. 

Numerical results. We start with a survey of the Nernst 
signal in the C4-symmetric case. v/T is shown in Fig. la 
as function of the chemical potential and the second- 
neighbor hopping t2, with t^ kept zero. While v/T is 
negative for all parameters, it is small for approximately 
circular-shaped Fermi surfaces, including the cases of 
small and large band filling for |t2/ii| < 1/2. (Note that 
v/T remains non-zero in the limit of zero filling, as both 
the lattice effects and the Fermi energy vanish in this 
limit.) In contrast, the Nernst signal is large near the 
line t2 = m/4 where the 2d van- Hove singularity is at the 
Fermi surface. 

We now consider the sensitivity of the Nernst sig- 
nal to a nematic distortion, i.e., a hopping anisotropy. 
Both Vxy/T and Vyx/T vary linearly with e, thus we plot 
the relative linear variation {l/v)d{vyx — Vxy)/de\^=Q in 
Fig. lb. As the Nernst signal itself, its variation is small 
near the band edges for 1^2/^1! < 1/2, while it is large 
near the van-Hove singularity. Remarkably, in a large 
regime of parameters the relative variation is 2. . . 5, i.e., 
the Nernst signal responds to anisotropies 2. . . 5 times 
stronger compared to the kinetic energy. The sign of the 
anisotropy, Vyx — Vxy, is robust and such that the (nega- 
tive) Nernst signal is enhanced if the thermal gradient is 
applied along the direction of stronger hopping. 

Hole-doped cuprates are located at —t^/ti « 0.2 . . . 0.4 
and ii/t\ = —0.5 . . .0, i.e., rather close to the van-Hove 
regime. We now discuss in more detail the thermoelectric 
response for such parameters. We choose ti — 0.38 eV 
and present results as function of hopping anisotropy 
for different (fixed) band fillings n and different i2,3- In 
YBa2Cu306+5, the bare hopping anisotropy extracted 
from band-structure calculations is |e| w 3 . . . 4% with 
ta < U (i.e. e < for a; along the crystalline a axis) - 
this effect is primarily due to the coupling between planes 
and chains.^** Nematic order is expected to strongly am- 
plify this anisotropy, and theoretical treatments have sug- 
gested anisotropies up to |e| = 20%.^° 

Fig. 2 shows the anisotropic elements of the response 
tensors a and d, together with the Nernst coefficient, for 
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FIG. 1: Evolution of the Nernst effect as function of chem- 
ical potential /i and second-neighbor hopping {—t2/ti). a) 
Nernst signal u/T in the isotropic case, in units of 2.45 x 
10^V/(K^T) X ro/s for ti = leV. b) Relative linear change of 
the Nernst signal with hopping anisotropy, d{i'yx — Vxy)/{i'de). 
In both panels, the white area corresponds to an empty band. 
For —tilt\ > 1/2 the band minimum moves to (0,7r), (7r,0). 
For ti — fi/4: the van-Hove singularity is located at the Fermi 
level, where both the Nernst signal and its sensitivity to 
anisotropy are maximum. (Very close to this line the data are 
inaccurate due to discretization errors.) By particle-hole sym- 
metry, the data for t2/ti > can be read off using /i — /i. 



t2/ti = -0.4, t3 = and band filling n = 0.86. While the 
anisotropy in the longitudinal conductivity (Jxx,yy is mod- 
erate, Fig. 2a, there is a substantial anisotropy in axx,yy. 
Fig. 2b. Remarkably, both anisotropics constructively 
interfere to generate a large Nernst anisotropy. Fig. 2c. 
The signal peaks at some finite e where the Fermi-surface 
topology changes as a function of the anisotropy - this 
point again corresponds to a (now anisotropic) van-Hove 
singularity. Note that Vxy can even change sign at some 
finite e, leading to a formal divergence of Vyx/vxy (Other 
additive contributions to e.g., from pairing fluctua- 
tions, may render this statement insignificant.) 

A few remarks are in order. First, the thermopower 
S — OLxxl <^xx comes out to be negative. Fig. 2b, also for 
e = 0. This is in contrast to the common wisdom that 
a single hole-like Fermi surface leads to a positive ther- 
mopower. Indeed, we obtain a positive (negative) ther- 
mopower in the limits of large (small) band filling, but 
the sign change in the [i-ti plane is not tied to the topo- 
logical change of the Fermi surface (i.e. the van-Hove sin- 
gularity), but instead occurs above half filling for t2 < 0. 
Second, the sign of the anisotropy in cr. Fig. 2a, is robust. 
This is physically transparent: tx > ty implies a larger 
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FIG. 2: Anisotropic response functions for nematic order, 
with hopping anisotropy e. The hopping parameters are 
ti = 0.38 eV, t2 = -OAti, ta = 0, the band filling is fixed 
to n = 0.86. a) Longitudinal conductivity a^x, Oyy and 
b) thermoelectric tensor axx/T, ayy/T; c) Nernst coefficient 
Vxy/T, Vyx/T. The units are a [3.19 x 10^°/(^ni) x ro/s] 
(using the c-axis lattice constant of YBa2Cu306+i), a/T 
[2.36 X 10^^V/(K2nm) x ro/s], v/T as in Fig. 1. 



velocity along x, and hence generically Ux 



Third, 



the sign of the anisotropy {axx — Oiyy)i Fig. 2b, changes 
in the /i-i2 plane exactly at the van-Hove singularity. 
Notably, in Eq. (3) this sign change is compensated by 
the similar sign change of the Hall conductivity Uxy , such 
that sign of the Nernst anisotropy is robust. 

The influence of model parameters is illustrated in 
Fig. 3. For finite 1^/12 < 0, the van-Hove singular- 
ity is outside the regime of e and n which is relevant 
for cuprates. (This conclusion may be changed by bi- 
layer splitting, see below and Ref. 10.) Fig. 3 also makes 
clear that the response of v/T to moderate anisotropics 
is strongly non- linear near a van-Hove point. 

Discussion. The large Nernst anisotropy for a ne- 
matic state prompts the question how this compares 
to other states with rotation symmetry breaking. For 
cuprates, the prime candidate is stripe order. Following 
Ref. 18, we have also calculated the Nernst anisotropy 
in various stripe states with real-space period 8 and 16. 
In most cases, the Nernst anisotropy for realistic pa- 
rameter values is moderate, i.e., less than a factor of 
two. Exceptions are states with extremely elongated 
Fermi pockets, which we only encountered for period-4 
charge-only stripes (and for which no evidence exists in 
YBa2Cu306-i-5). Hence, stripe order is unlikely to ex- 
plain the large Nernst anisotropy in YBa2Cu306+5-"^ 



4 



a) I b) 




0.1 0.2 0.3 0.1 0.2 0.3 



Anisotropy E Anisotropy e 

FIG. 3: Nernst coefficient v jT as function of ttie flopping 
anisotropy e, for parameters t\ — 0.38 eV and a) t2 — — 0.32fi , 
ts = -0.5t2, n = 0.875, and b) t2 = -0.2ti, = 0, n = 0.875. 
Ttie units are as in Fig. 1. 

We thus confirm the interpretation of the data in 
Ref. 3, attributing the large anisotropy for 80K < T < 
150 K to a nematic state (while stripe-like order may set 
in at lower temperatures). Indeed, using e = —0.2 from 
Ref. 10, we can semi-quantitatively match the experi- 
mental results for YBa2Cu306+5: We obtain a resistivity 
anisotropy ffb/ca = <^yy/<^xx ~ 1-5 . . .2.5 (Ref. 11) and a 
Nernst anisotropy Vb/va = v^y/vyx ~ 4 ... 10 (Ref. 3). 

Let us briefly discuss effects which were neglected in 
our calculation. A vertical dispersion from inter-layer 
coupling will smear the 2d van-Hove singularity, i.e., will 
remove the sharp peak in the Nernst signal as function 
of filling or distortion. However, as the Nernst signal is 
enhanced and very sensitive to anisotropy over a large 



range of parameters, this will not qualitatively change 
our conclusions. For YBa2Cu3 06+5, the main effect of 
inter-layer coupling is a bilayer splitting of the disper- 
sion, which will move one of the van-Hove singularities 
closer to the physical parameter regime, likely enhanc- 
ing the Nernst signal. Finally, at elevated temperatures, 
the scattering rate will be both anisotropic and energy- 
dependent. While it is known for cuprates that the antin- 
odal scattering rates become much larger than the nodal 
ones, our conclusions would only be modified if the scat- 
tering near the two inequivalent antinodes becomes dif- 
ferent by a factor of order 5 - this we consider unlikely. 

Conclusions. Our key result is a surprisingly large sen- 
sitivity of the Nernst anisotropy to symmetry-breaking 
Fermi-surface distortions, rendering the Nernst effect a 
unique tool to detect electron-nematic order. The ef- 
fect is particularly strong for non-parabolic band disper- 
sion, as occurs near van-Hove singularities and/or half- 
filling. We have provided a concrete theory for the mea- 
sured Nernst anisotropics in YBa2Cu306+5,^''^ pointing 
towards electron-nematic order in the pseudogap regime. 
We propose to search for Nernst anisotropics in other 
candidate compounds for such order, such as Sr3Ru2 07. 
Our results in Fig. 2 suggest that anisotropics of the ther- 
mopower are worth investigating as well. 

We are particularly grateful to L. Taillefer for sharing 
data prior to publication. We thank S. Kivelson, A. J. 
Millis, S. Sachdev, and L. Taillefer for discussions, and S. 
Sachdev for collaborations on related work. This research 
was supported by the DFG through SFB 608 (Koln) and 
the Research Unit FG 538. 



^ V. Hinkov et al, Science 319, 597 (2008). 

^ L. Taillefer, Keynote fecture, M2S-IX Conference, Tokyo 

(Sep 9, 2009). 
^ R. Daou et al, preprint arXiv:0909.4430. 
" R. A. Borzi et al, Science 315, 214 (2007). 
^ S. A. Kivelson, E. Fradkin, and V. J. Emery, Nature 393, 

550 (1998). 

® S. A. Kivelson et al, Rev. Mod. Pliys. 75, 1201 (2003). 

M. Vojta, Adv. Pliys. 58, 699 (2009). 
* H. Yamase and H. Kohno, J. Phys. Soc. Jpn. 69, 332 and 

2151 (2000). 

^ C. J. Halbotli and W. Metzner, Phys. Rev. Lett. 85, 5162 
(2000). 

° H. Yamase and W. Metzner, Phys. Rev. B 73, 214517 
(2006). 

^ Y. Ando et al, Phys. Rev. Lett. 88, 137005 (2002). 
^ Y. Kohsaka et al, Science 315, 1380 (2007). 
^ S. A. Grigera et al, Science 306, 1154 (2004). 
" K. Behnia, J. Phys.: Condens. Matter 21, 113101 (2009). 
^ Y. Wang, L. Li, and N. P. Ong, Phys. Rev. B 73, 024510 
(2006). 

^ O. Cyr-Choiniere et al, Nature 458, 743 (2009). 



" A. Hack! and S. Sachdev, Phys. Rev. B 79, 235124 (2009). 
A. Hackf, M. Vojta, and S. Sachdev, preprint 
arXiv:0908.1088. 

A probe which directly couples to the bulk quadrupoiar 
order parameter of an electron-nematic phase is usually 
not available. 
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^xy — ^yx 
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